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Abstract 



We study the spectra of certain integro-differential equations aris- 
ing in applications. Under some conditions on the kernel of the integral 
operator, we describe the non-real part of the spectrum. 

MSC Primary: 45K05, 35P20. Keywords: Gurtin-Pipkin equa- 
tion, spectrum, Denjoy-Wolff theorem. 

1 Introduction 

The following integro-differential equations arise in several areas of physics 
and applied mathematics, namely in heat transfer with finite propagation 
speed [5], systems with thermal memory [TT] , viscoelasticity problems pQ 
and acoustic waves in composite media [I]. 

(i) Gurtin-Pipkin equations of first order in time (GP1) 



u tt (x, t) = au xx — I k(t — s)u xx (x, s) ds, a > 0, x G (0, 7r), t > 0. (2) 



u tt (x, t) = u xx + eu txx - k(t- s)u xx (x, s) ds, x e (0, 7r), t > 0. (3) 
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u t (x,t) — / k(t — s)u xx (x, s) ds, xG(0,7r) 
Jo 

(ii) and of second order in time (GP2) 



t > 





(iii) Kelvin- Voigt equation (KV) 
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Here 



POO 

k(t) = / e- tT dfi{r) 
Jo 



is the Laplace transform of a positive measure dfi. We identify this measure 
with its distribution function //, so /i is increasing, continuous from the right, 
and the integral is interpreted as a Stieltjes integral [2]. We always assume 
that k is defined and integrable on (0, oo), that is 

^ < oo, (4) 

and that \i is supported on (do, oo) with some d > 0. 

Remark 1 // equation ([T]) can be differentiated with respect to t, then we 
obtain a special case of (j2J).' 

u tt = k(0)u xx - k{t- s)u xx ds. 
Jo 

with 

k = -^k. (5) 
dt v ; 

Remark 2 In the case k(t) = const = a 2 , equation t^j is in fact an inte- 
grated wave equation. Indeed, differentiation of (pQ) gives 

u tt = a 2 u xx . (6) 

If k{t) = a 2 e~ bt , then differentiation gives a damped wave equation 

u tt = a 2 u xx - bu t . (7) 

Let the initial conditions be «(•, 0) = £ for ([T]), and «(•, 0) = £, u t (-, 0) = r\ 
for © and ©• 

First we apply Fourier's method: we set <p n = sin nx and expand the 
solution and the initial data in a series in (p n 

00 oo oo 

u(x,t) = ^u n (t)v?„(x), f(ar) = ^^„(a;), rj(x) = ^r] n (p n (x). 
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The components u n satisfy ordinary integro-differential equations 
(i) GP1 

u n {t) = -n 2 / k(t - s)u n (s)ds, t > 0. (8) 



(ii) GP2 

u n (t) = —an 2 u n (t) + n 2 [ k(t — s)u n (s)ds, t > 0. (9) 

(iii) KV 

u n (t) = —n 2 u n (t) — en 2 u n + n 2 k(t — s)u n (s)ds, t > 0. (10) 



o 





We will denote the Laplace images by the capital letters. Applying the 
Laplace Transform to ([8]), and ffTOl . and using the initial conditions, we 
find 

(i) GP1 

zU n (z)-U = -n 2 K(z)U n (z) (11) 

or 

U M = ■ %t v ( 12 ) 

z + n z K{z) 

(ii) GP2 

z 2 U n (z) - z£„ -r] n = -an 2 U n (z) + n 2 K(z)U n (z) 

or 

U »W ~ z i + an2-n*K{zy [ ' 

(iii) KV 

z 2 U n {z) - z£„ - r) n = -n 2 U n {z) - en 2 [zU n - f n ] + n 2 K(z)U n {z), 

or 

TT ( \ + ~ ^ n £,n (1 A\ 

Un \ Z > = ~2~i 2~i 2 2FTT 14 

z^ + ezrr + rr — rriv (z) 

Denote the denominators in ( IT2|) . (TTBl) . and (IT3|) by F n (z), G n (2), and 
H n (z) respectively: 

F n (z) = z + n 2 fsT(z), G n (z) = z 2 + an 2 - n 2 K(z), 

and 

H n (z) = z 2 + ezn 2 + n 2 - n 2 K(z). 

Let and be the sets of zeros of -F n (z), G n (z), and H n (z) respec- 

tively. Set 

oo oo oo 

AgPI = [J-^n; A G P2 = [J G° , A KV = [J 
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Definition 3 The sets A GP i, A GP2 , and A KV are the spectra of equations 
([T]), (j2J) and ([3]) respectively. 



Study of the spectra is important for applications. See [I] and the references 
therein. 

Remark 4 Suppose that in the integro- differential equations (Tj^, (G|) and 
we replace the zero lower limit in the integrals by — oo. Then X is a point of 
the spectrum of F n or G n or H n if and only if there is a solution of the form 

u x (x,t) = e A VnO)- 

For such systems a semigroup approach to the equations is possible, (V. V. 
Vlasov, private communication). 

Remark 5 If © holds, then K(z) = k(0) - zK(z). 

We do not study here the regularity of solutions and consider the solutions 
as sequences {u n (t)}. Regularity of GP1 is studied in [8] and of GP2 in [12]. 
In [5], under the assumption that k(t) is twice continuously differentiable 
it was shown, in particular, that the solution u(x,t) of ([T]) is a continuous 
L 2 (0, T)-valued function. In [12], the conditions on the kernel k(t) and the 
initial data are found such that there exists a strong solution: for a 7 > 



L 2 (0,7r) + Ikxxt/, t)\\ 

L 2 (0,7r) 



dt < 00. 



We single out the important special case: the case of discrete measure fi 
with atoms at b\~ > of mass > 0, 

00 

^(*) = £7^' 0<6 1 <...^+oo, a k >0, (15) 



k=l 



Discrete measures arise in applications [4], where parameters a^, b^ are con- 
nected with auxiliary boundary value problems arising under averaging. 

2 Main Results 

In the general case the non-real part of the spectrum is described as follows. 
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Theorem 6 (i) For every n, each set F®, G^, or H® contains at most one 
point in the upper half-plane, and this point, if exists, belongs to the second 
quadrant. 

(ii) For n large enough, the set G° n contains a point z n such that 

z n = iy/a n + o{n) . (16) 

(in) If 

poo 

A= dfi(t) < oo, (17) 
Jo 

then for n large enough, the set F® contains a point z n such that 

z n = i\J~An + o(n). (18) 

Under the additional assumptions on K, we can find more precise asymp- 
totics of the non-real part of the spectrum. 

Theorem 7 Suppose that 

fi(t) = bt p + 0(t a ), t^oo, (19) 

where < a < p < 1. Then: 

(i) for n large enough, the set F° contains a zero z n of F n (z) such that 

^ e «r/(a-p) n 2/(2-p) ( X + / n 2(a-p)/(2-p) ) X ? (2Q) 

sin7rp / 

(ii) forn large enough, the set G° n contains a zero z n of G n (z) such that 

Zn = % J~an + J^a^-V^-Ml + o(l)). (21) 
2 sm7rp 

In the case of a discrete measure with finite number of atoms the next 
theorem was conjectured by V. V. Vlasov and N. Rautiyan (private commu- 
nication). 

Theorem 8 Let \i be a measure with compact support, that is 

k(t) = / e~' T d/i(r), < d < d < oo. (22) 



do 



Then the set A^v contains finite number of non-real points. 



In the case ([15]) it is not hard to study the real spectrum of the systems, 
see [7^\ and the discussion below. The simplest versions of theorems [6](iii) 
and [7] are contained in [J], [5], [7]. In [TU], a study of the spectrum shows the 
lack of controllability of the system. 



lr The assertions of Theorem 1 in |7] are correct only for n large enough. 
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3 Proof of the Main Results and discussion 



Our main tools are the Schwarz Lemma and the Denjoy- Wolff Theorem (see, 
for example, [9]). 

Schwarz 's Lemma. Let f be an analytic function which maps the upper 
half-plane C + into itself. Then the equation f(z) = z has at most one solu- 
tion w, and if such solution exists then \f'(w)\ < \, unless f is an elliptic 
fractional-linear transformation. 

Denjoy— Wolff Theorem. Let f be an analytic function which maps C + into 
itself, and suppose that f is not an elliptic fractional-linear transformation. 
Then there exists a unique point w G C + U {00} such that the iterates f* n 
converge to w uniformly on compact subsets ofC + , the angular limit f(w) = 
\im z ^, w f(z) exists and satisfies w = f(w). Moreover, the angular derivative 
f'(w) exists and satisfies \f'(w)\ < 1. 

Angular limit means that z is restricted to any angle e < arg(z— w) < n—e 
if w G R, or e < arg z < n — e ii w = 00. 
Angular derivative is the angular limit 

f'(w) = ]im(f(z)-f(w))/(z-w) 

z—tw 

if w G R; if w = 00 then it is defined by the angular limit 

— 1—r = lim f{z)/z. 

J [OO) z->oo 

The point w in this theorem is called the Denjoy-Wolff point of /. If w G 
R U {00} is such a point that the angular limit lim^^ f(z) = w and the 
angular derivative |/'(w)| < 1, then w is the Denjoy-Wolff point. 
Proof of Theorem 6. 

For the equations F n (z) = 0, G n (z) = 0, and H n (z) = 0, we will show 
that each of them has at most one solution in the upper half-plane. This 
solution belongs to the second quadrant. 

The Laplace transform of k(t) is 

PCO POO 

K(z) = / e~ zx / e- tx d[i(t)dx 
Jo Jo 




This is called the Cauchy transform of the measure dfi. Condition fll]) ensures 
that the integral defining K is absolutely and uniformly convergent on every 
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compact in the z-plane that does not intersect the negative ray. So K is 
analytic in the plane minus the negative ray. 
Moreover, 

Im K(z)lm z<0, z G C\K_. (23) 

We rewrite the equation F n (z) = as z = f(z) := —n 2 K(z), then / maps 
C + to C + and by Schwarz's Lemma has at most one fixed point in the upper 
half-plane. If Re z > 0, then 

f°° t + Re z 

Re + n 2 K(z)) = Re z + n 2 / rz-dfitt) > 0, 

Jo \ t + z \ 

so the solution must lie in the second quadrant. 

For G n (z) = 0, we first prove that there are no solutions in the first 
quadrant. Indeed z 2 + an 2 maps the first quadrant into C + , while n 2 K(z) 
has negative imaginary part in the first quadrant. 

To prove that G n (z) has at most one zero in C + we consider a branch of 
the square root which maps the lower half-plane onto the second quadrant. 
Then the equation is equivalent to 

z = f(z):=n^K(z)-a), (24) 

because all solutions are in the second quadrant. Function / maps C + into 
itself, and thus by Schwarz's Lemma can have at most one fixed point in C + . 

Similar argument applies to H n (z) = 0. There are no solution in the first 
quadrant. Indeed, if 

2 i 2,2 2 tsI \ 

z + ezn +n =nK{z) 

and z is in the first quadrant, then the LHS is in C + but the RHS is in C_. 
So the equation is equivalent to 

z = f(z) :=mf>{K{z)-ez-l), 

because all solutions are in the second quadrant. Function F maps C + into 
itself, and thus by Schwarz's Lemma can have at most one fixed point in C+. 
This completes the proof of part (i) of Theorem [6j 

To prove part (ii), we first notice that 

K(z) ->■ as z — re l9 , r -)■ oo, (25) 

uniformly with respect to 9 for \6\ < it — 5, for any given 5 > 0. This will be 
expressed by saying that K — > as z — > oo non-tangentially. To show this 
we use the following lemma. 
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Lemma 9 If 

| arg z\ < n — 5, (26) 

then 

\z + t\ x \z\ +t, t > 0. (27) 
Proof of the lemma. First, \z + 1\ < \z\ + t. Second, 

+ = \z\\l + t/z\ > \z\ cos S, 

and similarly 

\z + t\= t\l + z/t\ > tcos5. 
Thus \z + t\ > (l/2)(\z\ +t)cos5. This gives fl27] . 
Now 

|jf W | < o Mt) 



\z\ +t 



in the sector (j26"j) . and we obtain (|25|) . 
We rewrite the equation G n (z) = as 



z n = inyayl — K(z n )/a. 
As i^(-z) — > by Lemma 9, we obtain (ii). 

Now we prove (iii). Condition (|T7|) permits to obtain asymptotics of ii': 
K(z) = A/z + o{\z\~ l ), z^oo, (28) 
uniformly with respect to arg z in | argz| < ir — e. Indeed, 



\K(z)-A/z\ 



Now we rewrite F n (z n ) = as 

1 



- of MT*' Wt) = o(l * l ~ I) - 



z n = -n I A— + o(z L ) 

V z n 

which gives (iii). This completes the proof of Theorem 6. 

To find out whether a solution in C + exists for a given n, we consider the 
case of discrete measure \i. In the case that /i has a finite support, 

N 



our arguments are elementary. In this case, K(z) is a real rational function 
with iV poles, all of them on the real line. Then F n is a rational function 
of degree N + 1 because it has an additional pole at infinity. So it must 
have N + 1 zeros in the complex plane. On each interval Ik = (—bk+i, —bk), 
1 < k < N — 1 there is one zero by the Bolzano- Weierstrass Theorem. The 
remaining two zeros can lie either one in C + and one in C_ or both on 
some interval Ik, (which then contains three zeros), or both on the interval 
Io = (— b\, 0). One can give examples of each possibility. 

Examples. 

1. Suppose that the measure d\x has two atoms, that is 

k(t) = 1 (e- + e"*) . 

Then it easy to check that 

has 3 real zeros (and no non-real zeros). The additional two zeros are in 
(-1,0). 

2. Let 

k(t) = e - ' + 200e- 50 '. 
It is easy to check that the equation 

1 200 
K { x) + x = H — + x 

w x + 1 x + 50 

has 3 roots on the interval [—50, —1]. 

In the case of a measure with finitely many atoms, the question whether 
the Denjoy-Wolff point belongs to the real line can be solved in finitely many 
steps by using the criterion that all roots of a polynomial equation are real, 
see, for example [3]. 

In the general case of a discrete measure \i with atoms at bk, we denote 
Ik = (—bk+i,—bk),k > 1, and I = ( — 6 1? 0). If there is a solution w in the 
upper half-plane, it must be the Denjoy-Wolff point of f(z) = —n 2 K(z). If 
there is no solution in the upper half-plane, then the Denjoy-Wolff point w 
belongs to some interval Ik, and that —1 < f'{w) < 1. 

So theoretically we can find out whether there is a solution in the upper 
half-plane, by iterating /, starting from any point in C + , for example from 
the point z = i. The sequence Zk = f n {zk-\), k — 0,1, ... , must converge. 
If it converges to a point in C+, then this point is the unique solution in C+. 
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Convergence in this case is geometric. If Zk converges to a point on the real 
line then there is no solution in C + , but convergence in this case may be 
extremely slow, if \f'(w)\ = 1. 

Let w be the Denjoy-Wolff point of /. If w G C+, then w is the unique 
solution of F n (z) = in C+. The case w = oo is excluded because the 
angular limit of f(z) as z — > oo is 0. If w G R and w G Ik for some k, then 
f'(w) G [—1, 1] and this I& is the unique interval of the Ij which contains two 
additional real zeros of F n . 

We conclude that in the case w G C+, each interval Ik, h > 1 contains 
one solution of F n (z) = while Iq contains no solutions. 

The situation with GP2 and KV are similar: it has a solution in C + if 
and only if the Denjoy-Wolff of the function / = n<p(K — a) is in C + , and 
the Denjoy-Wolff point can in principle be found by iteration. 

Now we prove the theorem [71 
Lemma 10 Under the assumption / TJPj) we have 

K( z ) = -^z p - 1 + 0(z a - 1 ), Izl^oo, 
sin up 

uniformly with respect to argz in any angle \ &rgz\ < it — 5. Here the use the 
principal branch of z p ~ x which is positive on the positive ray. 



Proof. First, 



dt p 



^zP-\ 



Iq z + t shi7rp 

see, for example, (3 Probl. 28.22] or [131 Probl. 878]. 

So it is sufficient to prove our Lemma for the case that fi(t) = 0(t a ). We 
integrate by parts: 

R^ooJq t + Z R^oo\z + R J (t + z) 2 J 
In a sector | arg^| < tt — 5 this gives 

\K(z)\<C\im (j4%+ r^TT^ 
R^oo\\z\+R J t 2 +\z\ 2 



oo j.a 
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This proves the lemma. 

Now the proof of Parts (i) and (ii) follows the scheme of the proofs of 
Theorem El Parts (i) and (ii). 

Proof of Theorem [HI 

In view of Theorem 6, it is enough to prove that all roots of H n (z) are 
real for n large enough. We rewrite the equation H n (z) = as 

n 2 f{z) = z\ 

where 

f{z) = K{z)-ez-l. 

We know from Theorem 6 that all solutions belong to the second quadrant. 
Function / maps the upper half-plane into the lower half-plane. Let <fi be 
the branch of the square root that maps the lower half-plane onto the second 
quadrant. Now we rewrite our equation as n<j)(f(z)) = z. Function / is 
strictly decreasing on (— oo, — d) and f(x) ~ — ex as x — > — oo. So there is a 
point r < —d such that f(x) > on (— oo, r]. We have n<fr(f(x)) ~ —n^J—ex 
as x — oo, It follows that 

n<p(f(x)) > x for x < x , (29) 

where x < depends on n, K, e. 

Now suppose that n is so large that 

n(f)(f(r)) < r. (30) 

This will hold for n large enough because 4>{f{r)) < as we established above. 
Comparison of (1291) with (1301) shows that there must be a point r G (— oo, r) 
such that 

n<j){f{r Q )) = r and — n(/)(f(x))\ x=ro G (0, 1]. 

This shows that tq is an attracting fixed point of the function n<ft(f), and 
application of the Denjoy-Wolff theorem completes the proof. 

We finish with the following 

Question. Can one extend Theorem 8 to arbitrary measure d\i satisfying 
© ? 
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